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Abstract 

It is commonly accepted that in concentrated solutions or melts high-molecular weight polymers display random- 
walk conformational properties without long-range correlations between subsequent bonds. This absence of memory 
means, for instance, that the bond-bond correlation function, P(s), of two bonds separated by s monomers along the 
chain should exponentially decay with s. Presenting numerical results and theoretical arguments for both monodis- 
perse chains and self-assembled (essentially Flory size-distributed) equilibrium polymers we demonstrate that some 
long-range correlations remain due to self-interactions of the chains caused by the chain connectivity and the incom- 
pressibility of the melt. Suggesting a profound analogy with the well-known long-range velocity correlations in liquids 
we find, for instance, P(s) to decay algebraically as s~ 3//2 . Our study suggests a precise method for obtaining the 
statistical segment length b e in a computer experiment. 
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1. Introduction 

Following Flory 's "ideality hypothesis" (0), one 
expects a macromolecule of size N (numbers of 
monomers per chain) in a melt (disordered poly- 
meric dense phase) to follow Gaussian statistics 
(0; 0)- The official justification of this mean-field 
result is that density fluctuations are small, hence, 
negligible. One immediate consequence is of course 
that the mean-squared size of a chain segment of 
curvilinear length s = m — n scales as 



R 2 (s) = ({r m=n+s - r„) 2 ) n = b e 2 s, 



(1) 



at least if the two monomers n and m are sufficiently 
separated along the chain backbone and local corre- 
lations may be neglected (s 3> 1). Here, b e denotes 
the statistical segment length ( 3) and the position 
vector of monomer i. Obviously, the chain end-to- 
end distance is then R e = R(s = N) = b e yN. 
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Recently, this cornerstone of polymer physics has 
been challenged both theoretically Q HJIH) and nu- 
merically (0; Q) for three-dimensional melts and ul- 
trathin films Ja;0). The physical idea behind the pre- 
dicted long-range intrachain correlations is related 
to the "segmental correlation hole" (0), p*(s) w 
s/R(s) d , of a typical chain segment of length s in d 
dimensions. Due to the overall incompressibility of 
the melt this is thought to set an entropic penalty 
SU(s) w p*(s)/p (p being the total density) against 
bringing two adjacent chain segments together (0). 
The swelling of the chains caused by this repulsion 
is most pronounced for small curvilinear distances 
where SU(s) is important but vanishes for large s 
(see Ref. |p) for logarithmic corrections relevant for 
ultrathin films, d — > 2 + ) at variance to the behaviour 
of dilute polymer chains in good solvent (0). 

In d — 3 dimensions, the segmental correlation 
hole effect is weak, SU(s) oc l/y/s, and a standard 
one-loop perturbation calculation (@) can be read- 
ily performed following closely Edwards (0). For in- 
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stance, for the chain segment size this yields to lead- 
ing order 
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R 2 (s) k r ^24/^ 

2 — = —j= With Kr = — ^ — 

b e S V s Pb e 



(2) 



if intermediate distances larger than the monomer 
(1 <C s) and smaller than the end-to-end distance 
(s <C N) are probed. The latter restriction avoids 
additional universal physics due to chain end effects 
and polydispcrsity which may be calculated in prin- 
ciple |6J). In the present study we do not focus on 
these effects. Note that the l.h.s. should be rigor- 
ously zero for an ideal chain with statistical segment 
length b e while the r.h.s. is proportional to the per- 
turbation potential SU(s). 

In this paper, Eq. (2) is checked numerically by 
means of a well-studied lattice Monte Carlo ap- 
proach and used to extract a precise value for 
the statistical segment size b e for asymptotically 
long chains (cf. Figs. 1 and 2). We compare melts 
containing only monodisperse chains with systems 
of equilibrium polymers (EP) l)9t llCf) where the self- 
assembled chains (no closed loops being allowed) 
have an annealed size-distribution of (essentially) 
Flory type. Finally, a striking demonstration of the 
intrachain correlations is obtained by the power-law 
decay of the bond-bond correlation function (cf. 
Fig. 3). We begin our discussion by briefly describ- 
ing the numerical algorithms which allow to sample 
dense melts containing the large chain sizes needed 
(N > 1000) for a clear-cut test of the theory. 

2. Computational questions 



cal moves through phase space are considered) . EP, 
on the other hand, have a finite and constant scis- 
sion energy E attributed to each bond (independent 
of density, chain length and position of the bond) 
which has to be paid whenever the bond between two 
monomers is broken. (Apart from this finite scission 
energy for EP all our systems are perfectly ather- 
mal. We set ksT = 1 and all length scales will be 
given below in units of the lattice constant.) Stan- 
dard Metropolis Monte Carlo is used to break and 
recombine the chains. Branching and formation of 
closed rings are explicitly forbidden, such that only 
linear chains are present. 

The monodisperse systems have been equili- 
brated using a mix of local, slithering snake and 
double pivot moves (p). This allows to generate 
ensembles containing about 10 3 independent con- 
figurations with chain length up to W = 4096 
monomers. We have sampled EP systems with 
scission energies up to E — 15, the largest energy 
corresponding to a mean chain length (N) w 6011. 
For EP only local hopping moves are needed in 
order to sample independent configurations since 
the breaking and recombination of chains reduces 
the relaxation times dramatically, at least for large 
scission-recombination attempt frequencies. In fact, 
all EP systems presented here have been sampled 
within four months while the N = 4096 sample 
of monodisperse chains alone required about three 
years on the same XEON processor. EP are there- 
fore very interesting from the computational point 
of view and allow for an efficient test of predictions 
applying to both system classes. 



For both monodisperse and equilibrium polymer 
systems we compare data obtained with the three di- 
mensional bond-fluctuation model (BFM) (Q) where 
each monomer occupies the eight sites of a unit cell 
of a simple cubic lattice. In all cases presented here 
cubic periodic simulation boxes of linear size L = 
256 containing 2 20 monomers have been used. These 
large systems are needed to suppress finite size ef- 
fects, especially for EP (floT) - The monomer num- 
ber relates to a number density p = 0.5/8 where 
half of the lattice sites are occupied (volume frac- 
tion 0.5). It has been shown that for this density 
the BFM has a low compressibility which is com- 
parable to real experimental melts (|8|). The chain 
monomers are connected by (at most two saturated) 
bonds. In monodisperse systems adjacent monomers 
are permanently linked together (at least if only lo- 



3. Numerical results 

We now demonstrate numerically the long-range 
intramolecular correlations predicted by theory. 
Data for both monodisperse and equilibrium poly- 
mers (open symbols) are presented. As one expects, 
the size-distribution of polymers is found irrelevant 
for small segments, s -C (N). Note that averages are 
taken over all possible pairs of monomers (n, m = 
n + s) and, hence, the statistics deteriorates for 
large s. For clarity and to focus the discussion on 
internal correlations without chain end effects, data 
points for large s have been omitted for most of 
the indicated samples. For our flexible BFM chains 
already small curvilinear distances of about s w 10 
are well fitted by the theory, provided the total 
chain length is large, (N) 3> 10 3 . 
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Fig. 1. Curvilinear distance R 2 (s)/s for monodisperse chains 
and equilibrium polymers (EP) obtained by Monte Carlo 
simulation of the bond-fluctuation model (BFM) vs. curvi- 
linear distance s. We use log-linear coordinates to emphasize 
the power-law swelling over several orders of magnitude of 
s. The data points approach the asymptotic behaviour from 
below, i.e. the chains are swollen, as predicted by Eq. (2). 
Provided s (N), both system classes are identical. For 
larger s intricate chain end effects become relevant — es- 
pecially for the polydisperse EP as indicated for the scis- 
sion energy E = 11 — which are out of the scope of this 
study. The averages are taken for all possible monomer pairs 
(n, m = n + s). The statistics deteriorates, hence, for large s. 




Fig. 2. Replot of the segment size R(s), as suggested by 
Eq. (2), for monodisperse chains of length N = 2048 for 
different trial statistical segment lengths, as indicated. The 
bold line for the value b e pa 3.245 shows that the correction 
is indeed proportional to the correlation hole, This 
procedure is very sensitive to the value chosen and allows 
for a precise determination. The same b e allows to collapse 
data points from all systems (not shown) if s <C (TV) and 
100 < (AT). 

The size R(s) of curvilinear chain segments of 
length s, as defined in Eq. (1), is presented in 
Figs. 1 and 2. The first figure clearly shows that the 
chains are swollen. As predicted, the asymptotic 
Gaussian behaviour R 2 (s)/s — > b e 2 (dashed line) is 
approached from below and the deviation decays 



Fig. 3. The bond-bond correlation function P(s) for both 
system classes considered (at number density p = 0.5/8) 
as a function of the curvilinear distance s. Provided that 
1 <C s <g (N), all data sets collapse on the power-law 
slope (bold line) predicted by Eq. (4). The coefficient Kp 
indicated is consistent with a statistical segment length 
be « 3.245 ± 0.005. 



as SU(s) oc l/\/s (bold line). (The more intricate 
non-monotonous behaviour for large s for both sys- 
tem classes, especially for EP, is indicated for one 
example with scission energy E = 11.) 

Computationally, R(s) is an important quantity 
since it allows for the precise determination of the 
statistical segment length b e by extrapolation to 
asymptotically large segment sizes by means of the 
one parameter fit suggested by Eq. (2). The bold line 
indicated corresponds to b e — 3.245 which nicely 
fits the data over several decades in s. Note that 
a systematic underestimation of the true statistical 
segment length would be obtained by taking sim- 
ply the largest R 2 (s)/s w 3.23 2 value available, say, 
for monodisperse chains of length N = 2048. A bet- 
ter representation for performing the fit is to plot 
R 2 (s)/s vs. l/\/s m linear coordinates (not shown). 
This has the disadvantage that data points with s 3> 
10 are less visible. A much more precise representa- 
tion is given in Fig. 2 where we plot in logarithmic 
coordinates for different trial values of b e the l.h.s. 
vs. the r.h.s. of Eq. (2), i.e. the non-Gaussian devi- 
ations against the depth of the segmental correla- 
tion hole. Since both sides of the relation must be 
identical (neglecting higher order perturbation cor- 
rections) the fit value should yield a data collapse 
on the bisection line (bold line) . Note that our best 
value yields a collapse over more than an order of 
magnitude without curvature. 

An even more striking violation of Flory's ide- 
ality hypothesis may be obtained by comput- 
ing the bond-bond correlation function, P(s) = 
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(l ro=n+s • 1„)„ /I 2 - Here, 1, = r i+1 - r 4 denotes the 
bond vector between two adjacent monomers i and 

1 + 1 and l(p = 0.5/8) ss 2.634, the (mean-squared) 
bond length. It can be shown that P(s) must decay 
exponentially if no long-range memory effects were 
present due to the multiplicative loss of informa- 
tion transferred recursively along the chain Q) . The 
bond-bond correlation function is in fact the second 
derivative of R 2 (s) with respect to s 

P(s)l 2 = -^R 2 (s), (3) 

hence, it allows us to probe the non-Gaussian cor- 
rections directly without the trivial ideal contribu- 
tion l|6j). (We remember that the velocity correlation 
function of a fluid is similarly related to the second 
derivative with respect to time of the mean-square 
displacement.) This relation together with Eq. (2) 
immediately yields an algebraical decay 

P(s) = ^ with k p = KR (b e /l) 2 /8 (4) 

of the bond-bond correlation function for dense so- 
lutions and melts, rather than the exponential cut- 
off expected from Flory's hypothesis. 2 This predic- 
tion is perfectly confirmed by the data collapse pre- 
sented in Fig. 3 for both system classes. In principle, 
b e may also be measured from the power-law am- 
plitude, however, to less accuracy than the previous 
method (Fig. 2). 

4. Conclusion 

In this study we have challenged the famous ide- 
ality hypothesis suggested half a century ago by 
Flory. Testing recent theoretical predictions (|5j;|6j) by 
means of extensive lattice Monte Carlo simulations 
of linear polymer melts having either a quenched 
and monodisperse or an annealed size-distribution 
we provide strong evidence that long-range corre- 
lations exist in dense polymer melts which tend to 
swell the chains. Their most striking effect is the 
power-law asymptote for the bond-bond correlation 

2 Interestingly, the scaling of the bond-bond correlation 
function is related to the long-time power-law correlations of 
the velocity-correlation function found in dense fluids nearly 
fourty years ago il ill . Albeit the constraints which cause 
the self-interactions are different (momentum conservation 
for fluid particles, incompressibility for polymer chains) the 
weight with which these constraints increase the stiffness of 
the random walker is in both cases proportional to the return 
probability, l/s d / 2 , of a random walker in d dimensions. 



function suggesting a profound analogy with the 
well-known long-range velocity correlations in liq- 
uids jllh . We finally emphasize that our results rely 
on generic physics which should apply to all poly- 
mer melts containing long and preferentially flexible 
chains. 
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